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Abstract 

In many new physics scenarios, the particle content of the Standard Model is extended and the 
Higgs couplings are modified, sometimes without affecting single Higgs production. We analyse 
two models with additional quarks. In these models, we compute double Higgs production from 
gluon fusion exactly at leading-order, and present analytical results in the heavy-quark mass ap- 



D . proximation. The experimental bounds from precision electroweak measurements and from the 



measured rate of single Higgs production combine to give significant restrictions for the allowed 
deviation of the double Higgs production rate from the Standard Model prediction as well as on 
the branching ratio for the Higgs decay into photons. The two models analysed eventually present 
a similar Higgs phenomenology as the Standard Model. We connect this result to the magnitude 
of the dimension six operators contributing to the gluon-fusion Higgs production. 
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I. INTRODUCTION 



The search for the source of electroweak symmetry breaking has dominated particle the- 
orist's efforts for decades. Now that a particle with many of the right properties to be the 
Higgs boson of the Standard Model has been discovered , the efforts turn to under- 
standing the properties of this particle. In the Standard Model, the couplings of the Higgs 
boson to fermions, gauge bosons, and to itself are firm predictions of the model. In models 
with new physics, however, these couplings can be different. 

The dominant production mechanism for a Higgs boson is gluon fusion, which is sensitive 
to many types of n ew p hysics. The simplest possibility is for new heavy colored scalars 

m 

and/or fermions js iS] to contribute to Higgs production. However, since the observed Higgs 
candidate particle is produced at roughly the Standard Model rate, extensions of the Higgs 
sector beyond the Standard Model are extremely constrained. For example, a model with a 
sequential fourth generation of chiral fermions predicts large deviations in the Higgs rates 
16|-l20| and is excluded by the limits on Higgs production for any Higgs mass below around 



600 GeV 



21 



22|. The properties of these potential new colored particles are further limited 



by precision electroweak measurements. Models in which the Higgs boson is composite 
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34l |. along with models whi ch g enerate new higher dimension effective operators involving 



the Higgs boson and gluons 
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361], can also induce a single Higgs production rate different 



from that of the Standard Model. Untangling the source of possible deviations from the 
Standard Model by measuring the production and decay rates of the Higgs boson will be 
quite difficult in models where there are only small differences from the Standard Model 
predictions. 

In this paper, we examine the extent to which the gluon fusion production of two Higgs 
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bosons can have a rate very different from that predicted by the Standard Model 
given the restrictions from electroweak precision physics and from single Higgs production. 
The observation of double Higgs production via gluon fusion is important in order to measure 
the cubic self coupling of the Higgs boson 39|, |40|. In the Standard Model, the rate is small, 



although 
are large 



;he 0(al) radiative corrections are known in the infinite top quark mass limit and 
42|. For a 125 GeV Higgs particle, the most likely channel for HH exploration 
43j, where studies have estimated that the LHC at full energy will 



41 



is gg — 7- HH — )■ 6677 

be sensitive to this process with around 600 fb~^. Using jet substructure techniques, the 
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HH —7- hhW^W and HH — )■ bbr^r channels may be available with about 600 fb ^ 44 1 



and 1000 fb" 



40|. This is clearly not physics which will be done during the early phase of 



LHC operations, unless the rate is significantly larger than in the Standard Model 45|. 
Double Higgs production can further be studied through vector boson fusion, which is 



46|. Vector boson fusion production of two 



also sensitive to the three Higgs self coupling 
Higgs bosons can be affected by new operators involving the W and Z gauge bosons and the 
Higgs, but is not sensitive to the new colored particles which contribute to the gluon fusion 
process. Hence the two production mechanisms can provide complementary information. 

Double Higgs production from gluon fusion first occurs at one loop and is therefore 
potentially modified by the same new heavy colored particles which contribute to single 



Higgs production. However, as pointed out in Ref. [36|, single and double Higgs production 
are sensitive to different higher dimension effective operators and in principle, the single 
Higgs production rate could be Standard Model like, while the double Higgs production 
could be highly suppressed or enhanced. Here, we consider the effects of both heavy vector- 
like and chiral colored fermions on the single and double Higgs production rates, and the 
interplay between them. We will not consider models with extended Higgs sectors, or with 
higher dimension non-renormalizable operators. 

For single Higgs production, it is useful to analyze the effects of non-Standard Model 
colored particles using a low energy theorem 47|] . The theorem can be formulated using the 



background field method in terms of the traces of the mass matrices of colored objects, which 
eliminates the need to diagonalize complicated mass matrices 48| . The low energy theorem 
can be extended to double Higgs production, where new features arise 3^. In models with 
extended fermion sectors (for example, in little Higgs models 49|-|55| ) there are contributions 
to double Higgs production containing more than one fiavor of fermion 56|] . These diagrams 
contain axial couplings to the Higgs boson which are non-diagonal in the fermion states and 
we demonstrate how these effects can be included using a low energy theorem. Low energy 
theorems are extremely useful for single Higgs production and generally give estimates of 
the total cross section which are quite accurate. For double Higgs production, however, the 
low energy theorems provide an estimate of the total rate which typically disagrees with 
the exact rate by 50% or more. The low energy theorem does not reproduce kinematic 
distributions accurately, but instead predicts high energy tails which are not present in the 
full theory 
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FIG. 1: Feynman diagrams for gg — >• HH in the Standard Model. 



In this paper, we study the effects of heavy colored fermions on the gluon fusion double 
Higgs production rate and show that agreement with single Higgs production requires the 
double Higgs rate to be close to that of the Standard Model. We demonstrate how this 
can be understood in terms of the effective operator approach of Ref. [ssl and discuss the 
limitations of the low energy theorem for gg — )■ HH. Interestingly, composite Higgs models 
and little Higgs models receive potentially large corrections to the gg — )■ HH process from 
the non-renormalizable operator ttHH. The observation of such a large effect would be a 
"smoking gun" signal for such models 0, [s^ Q. 

II. DOUBLE HIGGS PRODUCTION 

1. The Standard Model 

In the Standard Model, double Higgs production from a gluon-gluon initial state arises 
from the Feynman diagrams shown in Fig. [H The result is sensitive to new colored objects 
(fermions or scalars) in the loops and to the Higgs trilinear self-coupling. The amplitude for 

g'''^ipi)g''''{p2) ^ H{p3)H{p^) is 



PnPuP2)Fi{s, t, M, rrit) + {puP2,p^)F2{s, t, u, ml\ 



where Pi and P2 are the orthogonal projectors onto the spin-0 and spin-2 states respectively. 



jM'^/^, ^\ _ _ P1P2 



Pr{Pi^P2) = 9^ 

Pi ■P2 

2 

Pr{Pl,P2,P3) = g^" + — {mHPiP2 - 2P1.P3P2P3 - 2p2-P3PlP3 + ^ P3P3) ' (2) 

spj^ 

s, t, and u are the partonic Mandelstam variables, 

s = {pi+ P2f , t={pi- pzf , u = {p2- pzf , (3) 
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Pt is the transverse momentum of the Higgs particle, 



p| = , (4) 
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(5) 



and V = (v^Gf)"^/^ = 246 GeV. The functions Fi and F2 are known analytically 
Finally, the partonic cross section is given by 

da{gg^HH) _ al |Fi(s, t, m?)p + |F2(s, t, u, m^)P 
It " 2i%3^ ^2 ' 

where we included the factor of | for identical particles in the final state. 
In the Standard Model, the chiral fermions are 



^1= , rffi, (6) 




where ? = 1, 2, 3 is a generation index and the Lagrangian describing the quark masses is 

- = E >^'^L^d^B. + >^Tl^u^r + h-c. . (7) 

i 

Here $ = is the Higgs doublet, <l = i(T2$* and 0° = Note that in the 

Standard Model the Higgs couplings A"'"^ are purely scalar. In the following we will focus 
on the third generation quarks and use the standard notation = t, = b, with A3 = Ai 
and A3 = A2. 

In the Standard Model, the dominant contributions come from top quark loops. Analytic 
expansion of the amplitudes in the limit » s yields the leading terms 



F^'Us.t.u.m;) = Vs< IH , 

Fl'^\s,t,u,ml) = --J1 + --I + ^ ^-^ ^ + — 



2^ 11 pit. Q2ml-5s ,^fs'^\} 
^ 45 mf { 154m? \mf J J 



The leading terms in the inverse top mass expansion of Ec 



theorem" result and give the mt-independent amplitudes |37|, 



(8) 

[8] are called the "low energy 



(4 4m?r \ 

F2(s,t,u,m?) \let -> . (9) 



From Eq. [HI we can clearly see that the triangle diagram has no angular dependence and 
only makes an s-wave contribution. This result is expected since the triangle diagram has a 
triple-scalar coupling, which has no angular momentum dependence. For the box diagrams, 
at the lowest order in F2°^ there is angular momentum dependence reflected in which 
is expected from the spin-2 initial state and spin-0 final state. At (9(m^^) in F^°^ there is 
also an angular momentum dependent piece proportional to p\. Since the initial and final 
states for the Fi contribution are both spin-0, this is a somewhat surprising result. To gain 
insight into the angular dependence of F^"^ and further insight into -Ff"^, the functions can 
be decomposed into Wigner rf-functions, d{. ,,^, where j is the total angular momentum and 
Si (sf) is in the initial (final) state spin: 



Fl'-(s,t,u,m?) = ~s 



7 m]j 540m|^ - 116m|^s + 59s" \ g 



13s^-52m|s , 
+ 5040m^ "^'^'^^^^^[m) 



11 s — Amir 



45 



m 



62m% — 5s / s^ 
1+ .."^ o +0 —r 



154m^ ' V fni 



4o(^) • (10) 



Here 6 is the angle between an initial state gluon and final state Higgs, 



t = ml --{1-/3 cose) and 13=^1 (11) 

In F^°^, we can see the expected spin-0 s-wave component, d^Q, and an additional spin-0 
d-wave component, c/qq, at 0{m^'^). The s-wave and (i-wave components are orthogonal. 
Hence any angular independent observables, such as total cross section and invariant mass 
distribution, are independent of the component of F^°^ up to C(m^^). Finally, F!^"^ is 
wholly dependent on the initial state spin-2 rf-wave function rf^o; as expected from Eq. [H 

In Fig. [21 we compare the total cross section for double Higgs production at different 
orders in the large mass expansion against the exact result^, as a function of the center 
of mass energy in pp collisions. We use the CTIO NLO PDF set ^\ and run the strong 
coupling constant through NLO from its value as{mz) = 0.118. We fix = 173 GeV and 
rrib = 4.6 GeV. The low energy theorem results are quite sensitive to the scale choice, and 
typically reproduce the exact results to within roughly 50% error. This "agreement" between 
the infinite mass approximation (LET) and the exact result is not improved by the inclusion 



^ The exact result always includes the contributions from both the top and bottom quarks. 
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FIG. 2: Double Higgs production cross section as a function of the hadronic center of mass energy 
-v/5 in the infinite top mass approximation, LET, (sohd hues) and retaining the 0(-^) corrections 
(dashed lines), normahzed to the exact result. The black (red) curves choose as the renormalization 
and factorization scales /x = 2mH (/^ = Mhh = \/^)- 



of higher orders in the large mass expansion. In single Higgs production, the reliability of the 



infinite mass approximation has been investigated through NNLO [59|-|62|. Because of the 
shape of the gluon parton luminosity, which peaks at large values of x = m?^/ s and decreases 
rapidly, the largest contribution to the hadronic single Higgs cross section comes from the 
region below the top quark threshold, s < 4m^, where the large top mass approximation 
holds. As a consequence, finite mass corrections to single Higgs production have an effect of 
less than 1%. On the other hand, for double Higgs production the partonic energy is always 
s > Amfj and the condition for validity of the low energy theorem, s ^ 4m^, is typically not 
satisfied. 

Fig. [3] shows the sensitivity of the results to the choice of the PDF sets. The exact result 
has a small sensitivity to the choice of LO vs NLO PDFs. However, the infinite mass limit 
(LET) of the result is quite sensitive to the choice of PDFs. Including higher order terms 
in the top mass expansion does not reduce this sensitivity to the choice of PDFs. 

The inadequacy of the infinite mass approximation for double Higgs production becomes 
even more apparent when looking at kinematic distributions js?]. Consider for example the 
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FIG. 3: Total cross sections for HH production using CTEQ6L LO PDFS and CTIO NLO PDFs. 
The renormalization/factorization scale is /i = 2mH in (a) and fi = Mhh = \^ in (b). For all 
curves, is evaluated at NLO. 



invariant mass of the HH system, 

da{pp HH) 2Mhh 



dM, 



HH 



s 



Hgg^HH)^ 
ar 



dC 



gg 



dr 



(12) 



where S is the hadronic center of mass energy squared, Mhh = y/s, and t = In Fig. H] 
we analyse the impact of the finite mass corrections to the invariant mass distribution at 
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FIG. 4: Invariant mass distributions for Higgs pair-production at \fS = 8 TeV and ^/S = 14 TeV, 
for terms in the large mass expansion up to 0{m^'^) (Eq. ED and with the fuU mass dependence. 
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FIG. 5: Transverse momentum distribution for double Higgs production cross section. The Stan- 
dard Model exact result, the LET and the heavy top mass approximations up to are 
shown. We choose as the renormalization and factorization scales /x = Mhh = \fs and use the 
CTIO NLO PDFs. 



the \fS = 8 TeV and = 14 TeV LHC The inclusion of the 0{m^'^) corrections does 
not significantly improve the low energy theorem results. The m^^ terms fail entirely in 
reproducing the exact distribution, in particular at large values of Mhh- Similar features are 
observed in the Pt spectrum shown in Fig. [51 Even for very small Pt ttli, the infinite mass 
spectrum does not reproduce the distribution accurately, although the transverse momentum 
distribution is well described when including the Oirri^^) terms. However, for px > rrit, the 
results from the heavy mass expansion drastically fail to approximate the exact distributions. 
A similar behaviour has been observed for the differential cross section da/dpT in higher 
order corrections to single Higgs production [63 1. 



2. Non-Standard Model bottom quark Yukawa coupling 

We briefly discuss the role of the bottom quark loops which are omitted when using the 
low energy theorems. In Fig. [6l we show the exact kinematic distribution for double Higgs 
production in the Standard Model, along with the result of the low energy theorem. The 
bottom quark contribution is negative but negligible in the Standard Model {Cb = ^ = 1 is 
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FIG. 6: Invariant mass distribution for Higgs pair production at ^/S = 8 TeV, in the infinite top 
mass approximation (solid black), with the full dependence on mt, but no b quark contribution, 
(red dotted) and including bottom-quark effects for increasing values of the Higgs-bottom quark 
Yukawa coupling (dashed lines). 



a rescaling factor of the bottom Yukawa coupling i/hh with respect to the Standard Model). 
The result of the destructive interference between the top and bottom quark loops remains 
small even when the bottom Yukawa is scaled up by a factor of 10. Only enhancements by 
factors as large as 50 cause the bottom loops to dominate and give significant deviations. In 
the Standard Model (with CTIO NLO PDFs and fi = Mhh), at = 8 TeV, the infinite 
mass approximation for the two Higgs cross section is about 70% of the exact two Higgs 
cross section. This remains roughly true if Cb is increased to ~ 10. However, if the b quark 
Yukawa coupling is increased by a factor of 50, this ratio goes to 9, and the low energy 
theorem is wildly inaccurate. 



3. Additional heavy quarks 



A simple extension of the Standard Model with additional quarks of charge | which can 



mix with the Standard Model like top occurs in 



49 



-52 



little Higgs 



heavy charge —\ quarks 



64| and composite Higgs 



23 



many 



-26 



new physics scenarios, for example 



65 



28|-|34| models. There can also be new 



and the formulae in this section apply to both cases. We 
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FIG. 7: Additional Feynman diagrams contributing to gg — )• HH in models with new heavy quarks 
coupling to the Higgs boson through non-diagonal Yukawa interactions. 

will take the new quarks to be in the fundamental representation of the color group. For an 
overview of the latest lower bounds on the masses of the additional quarks, see for example 
Refs. j^, Ig^] ■ Note however that the experimental analyses always assume the new quarks 
to decay entirely either through W or though Z. This is not the case in our models, and 
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68 



69|. 



the experimental limits are therefore weakened 

In addition to the diagrams of Fig. [H where any of the heavy quarks can be running in 
the loop, the double Higgs production receives contributions also from the mixed diagrams 
with two different quarks of Fig. [71 The mass terms and the interactions with a Higgs boson 
of a pair of mass eigenstate quarks /j, fj (of the same charge) are 



H 



with 



Y 



E7 



mi6ij + —H ) fj^R + h.c 



moij + — -H fj , 



A, 



Vij - Vji 



(13) 



(14) 



2 ' 2 

We consider real couplings. Therefore Yij = Yji and Aij = —Aji, and only the terms involving 
two different quarks fi and fj contain pseudo-scalar couplings, 

'Yij + -f^Aij 
Jiy- 

In the Standard Model Yu = rrii and Aij = 0. 
For arbitrary masses rrii and rrij, 



(15) 



F*"(., t, ml m]) = ^Fr\s, t, u, m^) + ^F^^s, t, u, m]) 

|-.IV(,A.-)+0(il). 



rrii 
Ami 



(16) 
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where y and M. are the Yukawa and the heavy quark mass matrices from Eq. 
box topologies, the leading terms in the large quark mass expansion are 



For the 



4 



+ 



rrr 



-\sTT{yM-'yM-')+o{£-^ 



F2{s,t,u,mi,mj) 



(17) 



The relative minus sign between the vector and axial contributions comes from Eq. [T4l 

Although the leading terms of the triangle and box diagrams were calculated in the 
diagonal mass basis, the cyclicity of the trace and the fact that both Ai and y rotate 
according to the same unitary transformations allow one to cast the results in Eqs. ITGlandfTTl 
into a basis independent form. Hence the Yukawa and mass matrices can be evaluated both 
in the mass basis, where Ai is diagonal, and in the current basis. In the current basis, 
y = The infinite mass limit of both the triangle and box diagrams can also be obtained 
via the low energy theorems 
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In our calculations in Sections IIII Al and llll Bt we retain the full dependence of the leading 
order amplitude on the quark masses. However, for small mass splitting 6 = m'j — mf the 
sub-leading terms have a simple and useful form. 



y2 _|_ y2 I 2y2 

" V '' F^^^is, t, u, m^) + ;-s^^ 
mf 6 mf 



4 Y^, + Y^^ S 



2 r 



8 Ai 
3 



1 + 



4s 



6 



GOmj 



2mj 



+ — ,— 



1 + 



7 m]j 
10 7^ 



mf 



mf 



22 pI 
45 m'f 



3 m^ m^ \m* / 



(18) 



Following 70|, we consider the infinite quark mass limit of these results and recast them 
into a convenient form for the calculation of the amplitudes for single and double Higgs 
production in models with extended quark sectors with respect to the Standard Model 
amplitudes. In the infinite mass approximation, the leading order amplitudes can be written 
as (Eqs. M E 



Agg^H oc Tr {yM-') 



A%%HH oc Tr {yM-'yM-') , 



(19) 
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where the omitted proportionahty terms do not depend on the masses and Higgs couphngs 
of the quarks. In the Standard Model, yu = rrit. The amphtudes only depend on the omitted 
proportionality factors, which therefore cancel when taking the ratio to the Standard Model 
result: 



gg^H 



■pbox 

^gg^HH 



^ = Tr{yM-') = ^ilogdetM) , 

A box 

A 



In Eq. [20] we used the relation y 



box,SM 
gg^HH 

dM 

dv 



Tr (yM-^yM-^) 



(20) 
(21) 



TOj. Eq. [2T]is equivalent to the result of Ref. [34 1 



III. EXAMPLES 



A. Singlet top partner 



We are interested in examining possible large effects in two Higgs production from gluon 
fusion in models which are consistent with precision electroweak measurements and the 
observed rate for single Higgs production. Topcolor models [23|, |28|], top condensate mod- 



els 24l-l27|. and little Higgs models 49l-l55| all contain a charge | partner of the top quark. 
We consider a general case with a vector SU{2)l singlet fermion, 7"^, which is allowed to 
mix with the Standard Model like top quark, 



68 



69|, 



7ll-l73| . The fermions are. 



R ■ 



(22) 



Following the notation of [5|, the mass eigenstates are t,T and b = (where t,b are the 
observed top and bottom quarks), and can be found by the rotations 



The chirality projectors are Pl,r 
parameterized as. 




^L,R 



' L,R 
'L,Rj 



(23) 



-J^-^ and the mixing matrices Uj^, f/jj are unitary and 



ui = 



cos 9l — sin 9 




- sin 9r 
cos 9fj 



(24) 
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We will abbreviate sl = sin 9l, cl = cos 9l- 
The fermion mass terms are 



where 



(1) 



=2 



A275 A375 



A. 



A, 



(25) 



(26) 



Without loss of generality, the T term can be rotated away through a redefinition of 
the right handed fields. The model therefore contains three independent parameters in 
the top sector, which we take to be mt,Mj- and 6l. The consistency of the model with 
electroweak precision measurements and its decoupling properties have been studied in many 



works 



m 



69 



7ll-l73| ■ We will not repeat this analysis here, but use the results of Ref. js]. 



It is interesting to note that in the limit 6*1, ~ (required by precision electroweak data), 
the mass terms for the top like quark and its partner become 



A, 



A, 



l + y(r-l) 



Mt 



1 + 



41 



where r 



Ml 



m 



. Decoupling of the heavy quark therefore requires s\ 



(27) 

as it was shown 



Since we are interested in Higgs production from the quark loops, we need the couplings 
to the Higgs boson, 

- Ch,1 = —CutLtRH + ^CttTlTrH + ^CtTtLTRH + "^CTtTLtRH + h.c. , (28) 
V V V V 

where 

Ctt = c\ , CtT = sl , CtT = CTt = SlCl ■ (29) 

Using Eq. [22] and the low energy theorems of Eqs. [20] and [2T1 it is straightforward to see 
that the single and double Higgs production rates are the same as the Standard Model 
up to corrections of O ^^j^grj- These corrections are further suppressed by the small 
mixing angles allowed by the bounds from electroweak precision data ^|. Both total and 
differential distributions are very close to the Standard Model (Fig. [8]), and one cannot use 
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FIG. 8: Invariant mass distribution in the Standard Model and in the top-singlet partner model 
(with Mt = 1 TeV) at the VS = 8 TeV LHC. 



double Higgs production to obtain information about additional vector singlet quarks. Fig. [H] 
uses the largest mixing angle allowed by precision electroweak data, and the reduction in the 
total cross section for the singlet top partner model from the exact Standard Model result 
is roughly 15%. This is of similar size to the reduction in the gg ^ H rate found in Ref. jsj. 
This model is an example of a case which will be extremely difficult to differentiate from 
the Standard Model. 



B. Mirror fermions 



As a second example, we consider a model which has a generation of heavy mirror 



fermions 



71 



74j-l77l|. There are four new quarks T\ and B^, B^, with charge | and 



— respectively. The quarks are in the SU{2)l representations. 




tLb\ 



R I'-'R 1 



R 



n 



q-2 »2 



(30) 



The first set of heavy quarks has the quantum numbers of the Standard Model quarks, while 
and B^ have have the left- and right- handed fermion assignments reversed from those 
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of the Standard Model. For simplicity, we assume there is no mixing between the heavy 
mirror fermions and the Standard Model fermions. This assumption eliminates the need to 
consider limits from Z ^ bb 65| and relaxes the restrictions from precision electroweak data 
discussed in Sec. IIII B iK 

The most general Lagrangian for the interactions of the mirror fermions with the Higgs 
doublet is, 



-C 



—2 



-7-2 ~ 



+XE^IJLrR + XfTrTI + XcBi^Bi + h.c. 



X5i + h.c.. 



(31) 



The mass eigenstates Xp {P = L, R; q = t,b) are obtained through unitary rotations 



cos 9p — sin 9p 
sin 9p cos 9p 



(32) 



and the mass matrices are 
Mu = 



Xb-^ 



Xf 



M 



D 



Xg Xc-^ 




(33) 
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We will denote the two top- like and the two bottom- like mass eigenstates as Ti,T2 and 
Bi, B2 respectively. The Lagrangian parameters Aj can be expressed in terms of the physical 
quark masses and the mixing angles. We report these relations in Appendix \^ 

Since all the quarks have different quantum numbers, it is not possible to rotate away 
any parameter in the Lagrangian. However, the SU{2) symmetry requires that 



M 



U,12 



M 



DA2 ■ 



(34) 



and therefore 



cos sin 9\ - Mr, cos 9^ sin 9p = Mb^ cos sin 9^ - Mb, cos 9^ sin 0^ . 



This relation can be written as 



sin 9t + 



sin 9t 



Mb,+Mb, 



sin ^1 + 



Mb, -Mb, 



(35) 



sin 9l , (36) 



^ We will not explore UV completions of this model that can mediate the decay of the mirror fermions 
through higher-dimensional operators and prevent the new quarks from becoming stable. 
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where eT = ef±ef\ 

The couphngs of the fermion mass eigenstates to the Higgs boson are 



£H "J^lJ^lrp rp TT I "-12-12^7-1 rp tt i "-11-12^7-1 rp TJ 1 "-12-11^7-1 rp tt i 



CT2T1 7f^ 



2v 



2v 



2v 



2v 



-DilUirH H iJ2L-D2R-n H 7: -DlL-D2/?-n H -D2L-Dli?-n + n.C. 



2t; 



2t; 



2t; 



2i; 



(37) 



where 



CTiTi 



MtJI + COS {2e^j) COS (2^^)] + Mr^ sin [20^) sin (2^^) 



2M5 



cos^ 0* - ~ (sin^ - sin^ 9' ) 



2Mti 

f'Ti cos {291) sin (20^ - Mt^ cos (20^) sin (2^j 

]Vf „ _ Mr , Mt + Mt 

iviT, _im^ sm{29l) - ^ sin(2^i) , 



1CT2T2, CT2T1, 



z z 

(c-TiTi, CT1T2) with ^ Mt2,91 



(38) 



Similar expressions hold in the bottom sector. 

The couplings to the electroweak gauge bosons that are needed for the computation of 
the Peskin-Takeuchi parameters (Sec. IIIIB ip are reported in the Appendix. 



1. Higgs production using low energy theorems in the mirror model 

For single Higgs production through top quark and mirror fermion loops, the low energy 
theorem of Eq. |20] yields 

where we introduce the fractional difference A of the single Higgs amplitude from that of 
the Standard Model. 

Both for simplicity and because one expects large corrections to the oblique parameters 
for a large mass splitting within each chiral doublet, we assume Mt^ = Mb^ = M and 
Mt2 = Mb2 = M(l + 6). In this limit. 



Agg^H = Ai;iH 1 1 + 4 - ^ [(2 + 5) sin 9t - 6 sin 0^ [(2 + 5) sin 9t + S sin 9",] | 



(40) 
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where we impose (see Eq. [36]) 

{2 + 5)smet+6smeX = {2 + 6) sinOt + 5 sinOl . (41) 

Given the recent observations at the LHC, we are interested in the case when 
Agg^H ~ Agg^jj. One simple way to recover this hmit is to have 

et-l, et-l^ (42) 

which for single production gives^ 

Agg^H - Al^^^ {l - ^l} = ^tt^H {l - ■ (43) 

To get the Standard Model result for gg ^ H further requires either 5 ~ or 0^ ~ 0^ ~ 0, 
where the constraint on the right-handed mixing angle in the top sector arises from Eq. HU 
The result of Eq. |l3] can be understood by inspecting the Yukawa couplings in the limit 

ct^^^ = -ct,t2 = -M6 cos\e'^) = -M6sm\2e'ji) , 

CT,T.= CT,T, = -^sin(2^U = ^sin(4^^). (44) 

Similar relations hold for the charge — | sector. Hence, for 5 ~ or 9*^^ ~ the diagonal 
Yukawa couplings go to zero and only the top quark, with its Standard Model Yukawa 
coupling, contributes to single Higgs production. The off-diagonal couplings of the mirror 
fermions to the Higgs boson are slightly less suppressed, and could induce deviations in the 
double Higgs rate from that of the Standard Model. 

From the low energy theorem of Eq. the box contributions to gg — )■ HH production 
(including top quark loops) can be estimated. 



pbo. ^ p.bo.,SM ^ ^^^^^ . 

AMI 4M2^ 4M|^ 4M|^ 2Mt,Mt2 2Mb,Mb2 
^ ^ al-al + al-al 1 (aj - alf + (aj - ajf 

^2 1 + 5 ^4 (5 + 1)2 ' y'^^) 



This relation holds for small 5. For 5 = 0, Eq. |36] requires sin6'*_ = sin0^, and 
Agg^H = Agg^^u {\ + Acos^ 9^_) . This result can be easily understood from the Yukawa couplings, 
CTiTi = CT2T2 — M cos^ 9'L and CT2T1 = —CT1T2 = ^sin(26''L). Also in this case, the gg ^ H rate is 
identical to the Standard Model rate for 9^ — ^. 
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where we defined 



02 = (2 + S) sin 9t , ce^ = S sin 6'5j_ 



^4 = (2 + (5) sin^: 



lb 



(46) 



For ^: 



|t,6 



f, Eq.SS] yields^ 



A. 



1 + 5 



cos^ 9i + 



2(1 + 5)2 




(47) 



Note that does not contribute in the infinite fermion mass hmit. The terms proportional 
to cos2(^^^) come from the contributions of the off-diagonal fermion-Higgs couplings. For 
this simple choice of parameters, the same term governs the deviations from the Standard 
Model both in single and double Higgs production. 

We are interested in determining how large a deviation from the Standard Model 
gg — )■ HH rate is possible with a minimal deviation in the gg ^ H rate. With the assump- 
tion of no mass splitting within the mirror doublets, there are five independent parameters: 
the mass scale M, which drops out in the heavy mass limit for the Higgs production rates, 
the mass splitting between families, 6, and three angles. Using Eq. HQl we replace one of the 
angles with the fractional deviation A of the gg ^ H amplitude from that of the Standard 
Model, 



We require this deviation to be within 10% and the mass splitting 6 between the two mirror 
families not to be too large (0 < 5 < 1), since we expect electroweak observables to put severe 
bounds on 6. Under these constraints, we perform a scan over 6,A,et and O^. The values 
of these parameters for which Eqs. Hi] and HH] yield real solutions for 9\^9^_ are represented 
by the blue dots in Fig. [9] The red diamonds represent regions where the difference Abox in 
the double Higgs amplitude from the box topology is larger than 15%. 

In the following, we fix = | in order to focus on a region with large Abox, and analyse 
how double Higgs production depends on 9^^ and 5 for a Standard Model gg ^ H amplitude, 
A = 0, and for ±10% deviations from it, A = ±0.1. This analysis is shown in Fig. [10] for a 
heavy mass scale M = 800 GeV. To qualitatively understand the features of these plots, one 

4 In the exact (5 = limit the result reads F^""" = pbox,SM _ ^^^^2 gb_ _^ g^os^ Ot] ■ 




(48) 
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=Mg =M, =Mjj =M(1+S), 0<5<1 




-1 -0.5 0.5 1 

+ 

FIG. 9: Allowed regions in the 9^_,9\_ parameter space where deviations, A, from the Standard 
Model gg ^ H amplitude are below 10% and the mirror fermion masses satisfy < < 1. The 
other two angles are fixed through Eqs. ST] and H8l The red diamonds denote regions where the 
gg —7- HH amplitude from the box topology deviates from the Standard Model by more than 15%. 



can consider the limit of small deviations from the Standard Model single Higgs amplitude 
and small family splitting 5, 

+ A [-1 + 5^ cos' eX + 0(5^)] + C(A2, . (49) 

For almost degenerate mirror fermions (5 ~ 0) and small deviations in single Higgs produc- 
tion from the Standard Model case, (which occurs when 9\_ = ±f ), the dominant term is 
^box ~ —A. When single Higgs production is suppressed, double Higgs production is always 
enhanced, while for a slightly enhanced Higgs single production rate, double production can 
also be suppressed. For A = and small 6, double Higgs production is also enhanced. In all 
cases, the minimal deviations from Standard Model double Higgs production occurs exactly 
at 9^ = while the maximum deviation is at 

Oi = arccos (2^) + ^f^A + (9(A') . (50) 
Finally, we note that the results of this section can be written in terms of an effective 



^box 



(5^cos^^ 



4 nb 



-5(1 - sin 
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FIG. 10: Deviations from the Standard Model box amplitude, F\°'^ , as a function of 6\_ for = 
M = 800 GeV and four different values of the fractional mass difference 5 of the two mirror famihes, 
for a 10% deviation in the single Higgs production amplitude (left plot) and for the same gg ^ H 
amplitude as in the Standard Model (right plot). The blank regions on the curves are not allowed 
for 9\,9^_ to be real. 



Lagrangian, which for 5 = is 

a 



C 



eff 



'1+4 cos' - - (1 - 4 cos' + 8 cos^ Q\)^ 
V 2v'^ 



(51) 



2. Bounds from electroweak precision data 



The new mirror quarks carry electroweak charges, and therefore contribute to the self 



energies of the electroweak gauge bosons [7^ |7^, l78| . A convenient way to parametrize these 



effects is through the Peskin-Takeuchi parameters 



79 



80| 



aASp 
aATp 
aAUp 



4,2 „2 



UzziMl) - Uzz{0) - n^-,(M|) - 



Cy/Sw 



n^z(M: 



Yiww{^) nzz(o) 



Mir 



As 



w 



Ml 
Iiww{Ml,) 



n 



WW 



(0) 



UzziMi) - UzziO) 



Ml 



^,z{Ml) 
Ml 



— s 



2 n,,(M|) 



W' 



Ml 



(52) 



where Uxyip^) denotes the transverse part of the vacuum polarization amplitude evaluated 



at momentum p' and c' 



w 



M'i 



1 — s^. The couplings of the mirror fermions to the 
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electroweak gauge bosons are reported in the Appendix. 



We use the fit to the electroweak precision data given in Ref . 8l| , 

5 = 0.03 ±0.10 , T = 0.05 ±0.12 , t/ = 0.03 ±0.10 
with correlation coefficients 



(53) 



( 1.0 0.89 -0.54^ 



\ 



I 



0.89 1.0 -0.83 
-0.549 -0.83 1.0 

The reference Higgs and top-quark masses are mH,rei = 126 GeV and nit^rei = 173 GeV. 
We use ruH = 125 GeV and so we need to account also for the Higgs contributions to the 
electroweak parameters. Up to terms of 0(M|/m|^), they read 



1 



H 



log 



m 



H 



m 



H,ref 



AT, 



H 



IGttc^ 



log 



m 



H 



m 



H,Tcf 



AU, 



H 



. 



The Ax^ is defined as 



Ax' = ^(X, - X,)(a'),r.i(X, - X,) , 



(54) 



(55) 



where Xj are the central values of the electroweak parameters from the fit in Eq. [53J Xj 
are the contributions to these parameters from the new mirror fermions and from the Higgs 
loops, and af- = aiPijaj , with a, being the errors given in Eq. ESI 

We consider the case of no mass splitting within the doublets, while the fractional mass 
difference between the two heavy families is parametrized by 6, and focus on the regions 
of parameter space where we expect the largest deviations with respect to the Standard 
Model gg HH amplitude, while the single Higgs rate remains very close to the Standard 
Model value. Following the discussion in the previous section, we therefore fix = |, 
A = {—0.1; 0; 0.1} and choose M = 800 GeV. In Fig. [Tl]we show the 95% allowed regions 
in the {sin 6'^, 6} parameter space for the three values of A (red bands), along with the 
regions where the box enhancement is larger than 15% (blue diamonds). The experimental 
bounds typically require 6 to be small. In this limit, the electroweak parameters assume 
simple expressions. 



ASf 
ATp 
AUf 



N, 



c 



A 



307r 
Nc 



2S 

— + 4:5 smel- A (I + 6 smO^,) ±0(5^ 
6 



A2 [2 + 6 {sme1 + 2) +0{6^)] , 



Nc_ 



A^ [2 + 6sme^^ + 0{6^)] 



(56) 
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95% CL Allowed Region From STU Fit 95% CL Allowed Region From STU Fit 

m„=125 GeV, A=-0.1, 6 '=It/2 m„=125 GeV. A=0, 9 '=11/2 




sin sin 8^ 

(a) (b) 



95% CL Allowed Region From STU Fit 

m^=125 GeV. A=0.1. 9 '=It/2 




sm 6^ 
(c) 

FIG. 11: Red bands: 95% confidence level allowed regions from the fit to electroweak data for single 
Higgs amplitudes which are suppressed (a) /enhanced (c) by 10% with respect to the Standard 
Model amplitude, or equal to the Standard Model amplitude (b), for 9'l = ^ and M = 800 GeV. 
Blue diamonds: parameter space regions which allow an enhancement of 15% or more to the double 
Higgs rate from the box topology. Such a large enhancement is not allowed by electroweak precision 
bounds in the case of A = 0.1 (c). 

where Nc = 3. For 5 — )■ 0, — t- = f and A — )■ (Eq. |43|) . However, a large increase in 
the double Higgs rate from the box topology can be obtained only for large values of 6. In 
particular, for A = 0.1 the electroweak precision observables do not allow the mass splitting 
to be large enough to obtain a significant enhancement, consistently with the results from 

Fig.[Tni 
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FIG. 12: Differential double Higgs production cross section in the Standard Model and in the 
mirror fermion model for OL = M = 800 GeV. The single Higgs production amplitude with 
respect to the Standard Model is suppressed by 10% (a), equal (b) or enhanced by 10% (c). We use 
CTIONLO PDFs and ^ = Mhh = ^fs- The curves labelled Low Energy Theorem use the infinite 
mass approximation to the rate. 



3. Phenomenology of the Mirror Fermion Model and — )• 77 

Once the parameters of the model are constrained to reproduce the Standard Model single 
Higgs amplitude to within ±10% and to be allowed by a fit to the precision electroweak 
data, there is very little freedom left to adjust parameters. The differential cross section 
for gg — )■ HH is shown for allowed parameters in Fig. [T2] and it is clear that this class of 
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models does not allow for a large enhancement of the HH production rate. The exact cross 
sections include both Standard Model t and b contributions, while the low energy theorem 
curves include the infinite mass limit of the heavy quark contribution. The largest allowed 
enhancement is found for A = —0.1 and in this case, the total cross section, pp — )■ HH is 
enhanced by ~ 17% over the Standard Model rate. 

The mirror fermions also contribute to the rate for H — )■ 77^. We again consider each 
mirror family to be degenerate between the charge | and charge — | quarks, and the two 
families to be split by a mass difference M6. In the limit tuh << rrit, Mw, M 87|, 



m ^ 77) 



TiH ^ 77) 



SM 



1-15 

47 



CTiTi CT2T2 I 1 CBiBi I Cb^Bi 

2Mt, 2Mt2 4:\2Mb, 2Mb2 



1 - — [5 + sin^^ (Ssin^^ -Ssin^i)] 
-^5sin^^ (sin^^ - sin^i) + 0{5^] 



(57) 



where we impose only the angle relation from Eq. |4T] and expand for small 5. In the limit 
5 = (and therefore = 6*_ from Eq. Hi]), the branching ratio into photons cannot be 
larger than in the Standard Model. 

We relate the deviations in the photon decay branching ratio to the deviation A from the 
Standard Model single Higgs production rate^. 



T{H ^ 77) 



T{H ^ 77) 



SM 



47 

-1a 

47 
A^ 



A{5 + 1] 



{5 + 2) sin 61 + 5 sin d\ 
{{5 + 2) sin^*_ - 5 sin OXf ~ (6 + 2) sinOl -6sinel_ 
3(5 + 2) sin 01 12(5+1) 



1 + 



{6 + 2) sin et- 6 sin el 
6 + 1 



((5 + 2) sinOl -5sin0^)2_ 

(58) 



47 ((5 + 2) sin 0L -5 sin 0^)2 

Imposing only the bounds from electroweak precision observables, and performing a gen- 
eral scan over the input parameters 6,9'^,9'^,9t^ (fixing 9i_ through Eq. HH M = 800 GeV 
and 5 in the range {—0.5; 2}), we find that the Higgs branching ratio into photons can 
have large differences from the Standard Model predictions, with suppressions as large as 
90% and enhancements up to 10%. Requiring also the single Higgs production rate to be 



^ We c onside r only the contributions of heavy mirror quarks. Heavy leptons can also affect the H 
rate 



82 



mi- 



ll 



This result holds for arbitrary values of the parameters. 
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FIG. 13: Ratio of T{H — )■ 77) to the Standard Model Higgs branching ratio into photons for 



the points of Fig. 11(b), where A = and 0^ = §. We fix 5 = ±0.2, which is allowed from 
the electroweak fit for all the values of Larger deviations from the Standard Model — )• 77 
branching ratios arise outside this range of 6, in the regions where | sin^^Ll is close to 1. 



close to the Standard Model value puts severe constraints on these deviations. For a single 
Higgs production amplitude equal to the Standard Model prediction, the maximum devia- 
tion in the Higgs branching ratio into photons is ±5%. For the regions of parameter space 
of Fig. 11(a), where 6'^ = | and a —10% deviation from the Standard Model prediction for 



the gg ^ H rate is allowed, only small enhancements (up to +10%) of the if — ?■ 77 rate 
are allowed. For a +10% enhancement in the single Higgs rate over the Standard Model 



prediction (Fig. 11(c)), the branching ratio into photons deviates from the Standard Model 



prediction by at most by a few percent. We show how these deviations depend on the free 
input parameters 6, sin 9^ in Fig. [T3|, where we focus on A = and pick two values of 6 



which are allowed by the electroweak fit over all the range of 6^ (with = |). The clear 
conclusion is that the restrictions from precision electroweak data, combined with a single 
Higgs production rate close to the Standard Model prediction, do not allow for significant 
deviations in the if — )■ 77 rate in this class of models. 
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IV. CONNECTION TO GLUON-HIGGS DIMENSION SIX OPERATORS 



An interesting idea |36| is to combine single and double Higgs production to gain in- 
sights on the mechanism giving mass to the particles that contribute to these loop-mediated 
processes. Including contributions up to dimension-6 operators, the effective Lagrangian 
responsible for the Higgs-gluon interactions can be written as 

C = CiCi + C2O2 . (59) 

Particles whose mass arises entirely from renormalizable Higgs couplings induce an operator 

O2 = —G'^G'^^'''' log ( ^] ~ ^GIG'^'^'' (-- ^] . (60) 

If the particle receives contributions to its mass from other sources as well, an additional 
operator 

. . ^^Gl^G-'" (f + g) (61) 

arises. In the Standard Model cf^'^ = 0, cf*''^ = 1. The two operators contribute differently 
to Higgs single and pair production and the different rates in these channels constrain the 



coefficients ci and C2. Following 36|], one can derive these two coefficients in a background 



field approach. The Higgs field is treated as a background field, and the masses of the 



heavy partic^ 



theories 



47 



es become thresholds in the running of a^. Matching the low- and high-energy 



where M.{H) is the Higgs-dependent mass matrix and 56/ = 2/3 for fermions in the funda- 
mental representation of the color group. This yields the effective Lagrangian 



^'^^ = I^^^-^"'"'^ ^ • (63) 

We write the determinant of the mass matrix as 

det M{H) = [l + Fi{H/v)] X P(Ai,mi,w) , (64) 

where P is a polynomial of the Yukawa couplings Aj and fermionic masses rrii and in general 
Fi{H/v) = F{H/v,Xi,mi,v). If Fi{H/v) is such that i^'(O) = 1 + Fi{0), and all the higher 
order derivatives vanish before electroweak symmetry breaking, then the Higgs production 
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rates via gluon fusion in the heavy quark hmit are exactly as in the Standard ModeF. This 
is the case in the singlet top partner model, where Fi{H/v) = H/v and therefore ci,2 = cff^. 

Interestingly, one can have the same single Higgs production rate as in the Standard 
Model, but a different double Higgs rate, only for F"{0) 7^ 0. If also the first condition, 
F-'(O) = 1 + Fi{0), is not met, then the single Higgs rate is not Standard Model like. In 
such a case, we note that for Fi independent of Yukawa couplings and fermionic masses. 



the Higgs rates do not depend on the details of the fermion sector [3J] and deviations can 
arise only from changes to the Higgs potential. If Fi depends on the Yukawa couplings and 
fermionic masses, the Higgs rates will in general be related to these parameters. Such a 
situation occurs for example in the mirror fermion model. In this case 



(l-/3t)2 



l + TT 



(l-/3t)2 



We define 



eaf 



In terms of the physical parameters, 

4(1 + 5) 



(l-/3i,)2 ' 
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aIw72 



(65) 
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Q 



t,b . 



(66) 



(67) 



(2 + 5)2 C0S2 0'L - 52 C0S2 01 ' 

For /3b — )■ 0, c5 and go to twice the Standard Model value. In this limit, the vector 
contributions to the fermion mass matrix vanish, and the masses come entirely from elec- 
troweak symmetry breaking. Since there are two quarks, an extra factor of two arises. In C2 
one clearly sees the +1 contribution coming from the Standard Model top quark. 

The coefficients governing single and double Higgs production are then 

1 1 



ch 

ChH = Ci - C2 



Ci + C2 = 1 + 2 

-1 - 



+ 



1 - A 1 - A 
1 + A , 



l + (3b 



(68) 



;i-A)^ (1-/3.) 

The two rates depend on the two independent parameters Pt, Pb from the top and bottom 
sectors. Even if we require the single Higgs rate, gg H, to be close to the Standard Model 



^ For the purpose of this discussion, we only need -F'/'(0) = 0. Nonvanishing derivatives at higher orders 
only affect gluon fusion production of three or more Higgs bosons. 
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FIG. 14: The shaded red regions correspond to amplitudes for gg ^ H within ±10% of the 
Standard Model rate. 



value, 



ch = ci + C2^ cf{l + A) = 1 + A 



(69) 



we are left with an independent parameter that can yield completely independent variations 
in the double Higgs rate. 

In Fig. [m we show the regions of and Pt which reproduce the Standard Model Higgs 
amplitude to within A = ±10%. Imposing the constraint of Eq. |69]on the single Higgs rate 
in general constrains the double Higgs rate, gg — )■ HH, 



CHH -> 2ci - (1 ± A) . 



(70) 



In the singlet case, ci = and deviations in single and double Higgs rates must be of the 
same order of magnitude. In the mirror case, Ci can deviate from zero, which removes the 
close relationship between single and double Higgs production. 
In terms of the parameters of the mirror fermion model, 

8 5 - A 




1-A 



^A±A^ 



(2 ± 5 f cos^ el - 52 cos2 
l + ~5 



±0(A) 



(71) 



The term in the curly brackets correctly reproduces 1 ± Afeo^. from Eq. HUfor A = 0, ^* 



2 • 



A large effect in the double Higgs rate requires large Ci, and, in turn, /S^ ~ 1. This is seen 
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FIG. 15: Enhancement of the box contribution to gg — )• HH for a single Higgs amplitude within 
±10% of the Standard Model prediction. 

in Fig. [151 where we fix Ph to reproduce the single Higgs rate within 10% the Standard 
Model value. However, from Eq. |67] /3( — ?■ 1 implies 5— lor5— t-oo. These are not 
viable solutions. The first one corresponds to massless quarks. The second one requires 
non-perturbative interactions with the Higgs (large A^, Ad) for heavy quarks (large A^, Ai;'), 
as in Eq. [661 In the mirror fermion model discussed in this paper, large deviations in the 
gg — )■ HH rate do not occur. 



V. CONCLUSIONS 



We analysed double Higgs production from gg — t- HH in the Standard Model and in mod- 
els with additional heavy vector or chiral quarks. In the Standard Model, we compared the 
approximate results in the large top mass expansion with the exact cross section, and anal- 
ysed the dependence of the production rate on the choice of the renormalization/factorization 
scale /X and on the PDF sets. As is well known 42], 5^ , the low energy theorems fail to accu- 
rately reproduce both the total and differential double Higgs cross sections. The differential 
distributions are poorly estimated by the low energy theorems and predict a large tail at 
high invariant masses. The discrepancy is smallest for the scale choice fi = 2mH, yielding a 
10 — 25% difference from the exact calculation of the total rate. Further, the predictions of 
the large top mass expansion depend sensitively on the choice of PDFs. Inclusion of higher 
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order terms in the large mass expansion does not improve the convergence towards the exact 
results. 

We discussed how the combination of single and double Higgs production from gluon 
fusion might give insight into the mechanism giving mass to quarks. The parameters of 
models with new heavy fermions are strongly constrained both by the observed rate for 
gg ^ H and by precision electroweak measurements. In the case of a new heavy vector 
singlet quark, electroweak precision observables strongly constrain its mixing with the top 
quark |5[. The singlet needs almost to decouple from the Standard Model particles, and 
therefore deviations from the Standard Model in both the single and double Higgs rates are 
small. 

The situation is more interesting in the case of heavy mirror quarks which are not allowed 
to mix with the Standard Model fermions. The bounds from electroweak precision data 
still allow for the single Higgs production cross section to differ from the Standard Model 
predictions. However, after restricting the deviations in the gg H rate from the Standard 
Model rate to be small, the resulting double Higgs cross section and distributions become 
close to those of the Standard Model. The reason for this behaviour becomes clear in terms 
of the two dimension six operators Oi and Once we fix the single Higgs rate to be close to 
that of the Standard Model, large deviations in the double Higgs rate occur only if one of the 
mirror family becomes very heavy, with non-perturbative Higgs interactions, or very light, 
outside the range ttlh < 2mq where the operator expansion applies. In the mirror fermion 
model we also investigated the effects of the additional quarks on the Higgs branching ratio 
to photons. After the constraints from the observed single Higgs cross section and precision 
electroweak measurements are taken into account, the branching ratio — > 77 is always 
within 10% of the Standard Model rate. 

Therefore, in the two example of models with new heavy fermions which we studied, the 
constraints from the observed gg H rate, combined with precision electroweak data, do 
not allow large deviations of the gg — ?■ HH rate from the Standard Model prediction. 
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Appendix A: Electroweak parameters in the mirror fermion model 



We present here some useful formulae for the mirror fermion model. 
The parameters Aj appearing in the mass Lagrangian can be expressed in terms of the 
physical masses and mixing angles as 

Ab = cos el COS 6*^ + sin sin , 
v2 

A £) —j= = Mt^ sin sin cosOt cos 6*0 , 

Ag = Mt^ sin 6*^ cos 6^^ - Mt^ cos $1 sin , 

= Mt^ cos 6*1 sin 6*^ - Mt^ cos 6*^ sin 9^. (Al) 

Similar relations hold for the corresponding parameters in the bottom sector, with My. 
Mb, and O'p ^ O^p. 

The charged current interactions among quarks of charge Q and [Q — 1) are 

= ^ E {^"q^" [^^t^^ + ^^Pn] g^Q.i) } + h.c. , (A2) 



with 



1/^^^ = cos 61^ cos 6'^ = sin6'^cos6'[ 

^T^B, = cos sin V^^p^ = sin sin 

^r'Jsi = sin sine*, V^^p^ = - cos sin 

V^fsi = -sin^^cos^*, V^^p^ = cos e^ cos e*,. 



(A3) 



We can rewrite these relations as 

The neutral current interactions among quarks of charge Q are 



where 



^TiTi = cos^ ei X^^T^ = X^^Ti = sin ei COS A|; = sin^ ei 

X^^^^ = sin^e'j, X^^^^ = X^^j,^ = - sine*, cos X#^^^ = cos^ e*, 



(A6) 
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The same relations, up to an overall minus sign, hold in the bottom sector. In more compact 
form we can write 



^S = ±K)„K)^.. X5 = ±(f/«)_^(f/^')^^. (A7) 
where the plus sign holds in the top sector and the minus in the bottom sector. 
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